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Abstract 



we calculate the new physics contributions to seven measured decays B — * 
ir + 7r~ , Kir and Krj 1 in the general two-Higgs-doublet model (Model III). 
Within the considered parameter space, we find that: (a) the CLEO/BaBar 
measurement of B —* 7r + 7r~ decay prefers a small i^^O): F^ n (0) = 
0.25 ± 0.03; (b) the new physics enhancements to the penguin-dominated 
B —* Kit and B — > Kr]' decays are significant in size, ~ (40 — 70)% w.r.t 
the standard model predictions; and (c) the new physics enhancements can 
boost the branching ratios B{B — ► K + rj') and B(B — > K°t]') to be consistent 
with the data within one standard deviation, and hence lead to a simple and 
plausible new physics interpretation for the rj'K puzzle. 
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One of the main objectives of B experiments is to probe for possible effects of new physics 
beyond the standard model (SM). Precision measurements of B meson system can provide 
an insight into very high energy scales via the indirect loop effects of new physics. The B 
system therefore offers a complementary probe to the search for new physics at other hadron 
colliders JTJ,§. 

Up to now, CLEO, BaBar and Belle Collaboration 0-|6| have observed fourteen B u ^ 
meson two-body charmless hadronic decay modes 

B -> 7T ± 7T T , KtT, Krf, p ± 7T T , p°7T ± , WTT* K*T], ffV, QK* . (1) 

The experimental measures are generally consistent with the theoretical predictions based 
on the effective Hamiltonian with generalized factorization |7|||, with an exception of the 
so-called rj'K puzzle: the B — > Kr]' decay rates are much larger than that expected in the 
SM 0. 

The unexpectedly large branching ratios of B — > Kr]' was firstly reported in 1997 by 
CLEO §, and confirmed very recently by CLEO and BaBar B(B -> K+r/') = (75 ± 

10) x 1(T 6 ( average of the CLEO and BaBar result), and B(B -> K°T]') = (S9±i|±9) x 10" 6 ( 
CLEO ). The Krj' signal is large and stable, and clearly much larger than the SM predictions 
B[B — > Kr]') = (20 — 50) x 10~ 6 as given in Refs. 0,^,0. In order to accommodate the 
data, one may need an additional contribution unique to the rj' meson in the framework of 
the SM, or enhancements from new physics models beyond the SM to explain the B — > Kr]' 
puzzle jT2 . 



In a previous paper [UJ1, we considered the second possibility and calculated the new 



physics effects on the branching ratios of seventy six B — > hih 2 decay modes in the general 



two-Higgs-doublet models (2HDM's) ||13|Jl4| , and found that the new physics enhancement 
to the penguin-dominated decay modes can be significant. In this letter, we focus on seven 
B — > PP decays (where P refers to the light pseudo-scalar mesons) whose branching ratios 
have already been measured. We firstly find the constraint on the form factor F B7r (0) from 
the measured B — > 7r ± 7r T decay rate, and then check the consistency between the theoretical 
predictions and the data for the four B — > Ktt decay modes after including the new physics 
contributions in model III: the third type of two-Higgs-doublet models |I3| , [I4| . We finally 
calculate the new physics enhancements to the B — > Kr]' decays and study the effects of 
major uncertainties. 

On the theory side, one usually uses the low-energy effective Hamiltonian with gener- 
alized factorization |T^|T^J^,g] to calculate the two-body charmless B meson decays. For 
the inclusive three-body decays b — > sqq with q 6 {u, d, s} the effective Hamiltonian can be 
written as M, 



Gp I , 

72 



H eff (AB = l) = r± |£<~ (y ub v* s Q] + v cb v; s Qi) - v tb v* 



10 

X! CjQj + CgQg 
i=3 



(2) 



where the operator basis contains the current-current operators Qi,2, the QCD penguin 
operators Q3-6, the electroweak penguin operators Q9-10 and the chromo- magnetic dipole 
operator Q g , the explicit expressions can be found easily for example in Ref. 0. For b — > dqq 
decays, one simply makes the replacement s — > d. 

Following Ref. [0], we also neglect the effects of the electromagnetic penguin operator 
Q7 7 , and do not consider the effect of the weak annihilation and exchange diagrams. The 
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coefficients Cj in Eq.(|2|) are the well-known Wilson coefficients. Within the SM and at scale 
M w , the Wilson coefficients Ci(M w ) , ■ ■ ■ , Ci (M w ) at next-to-leading logarithmic order 
NLO) and C g {M\y) at leading logarithmic order (LO) have been given for example in Ref. 

0- 



In a recent paper |T7|], Chao et al. studied the decay b — > 57 in model III by assuming 
that only the couplings X u = \\tt\e l9t and A&& = |Af>b|e i6 ' b are non-zero 1 . They found that 
the constraint on imposed by the CLEO data of b — > sj can be greatly relaxed by 

considering the phase effects of \ u and A^. From the studies of Refs. PJPTH, we know that 
for model III the parameter space 



A« = 0, for ij 7^ tt, or bb, 



|A«| = 0.3, |A 



bi>\ 



35, # = (0°-30°), M H + = (200 ± 100)GeV, 



(3) 



are allowed by the available data, where 9 = 6bb — Ott- In this letter, we calculate the new 
physics contributions to seven B meson decay modes in the Chao- Cheung- Keung (CCK) 
scenario of model III W7\. Since the new physics corrections on the branching ratios of 
two-body charmless hadronic B u ^ decays in models I and II are small in magnitude |L1 1 , we 
do not consider the cases of models I and II in this letter. 

Following the same procedure as in the SM, it is straightforward to calculate the new 
7-, Z°- and gluonic penguin diagrams induced by the exchanges of charged-Higgs bosons 
appeared in model III [fTT]| . After taking into account the new physics (NP) contributions, 
the Wilson coefficients Ci(My/) i — 1, • • ■ , 10 at the NLO level and C g at the LO level can 
be written as 



C x {Mw) 
C 2 (M W ) 
C 3 (M W ) 

C 4 (M W ) 

C 5 (M W ) 

C 6 (M W ) 

C 7 (M W ) 
C 8 (M W ) 



11 a s (Mw) 35 a e 



1 

6 An 
11 a s (M w ) 

2 



18 4vr 



47T 

24tt 



E (x t )+E^ 



NP 



6n sin 2 6 W 

a s (M w 



2B (x t ) + C (x t ) + C, 



NP 



8tt 
a s {M w 

24tt 
a s {M w ) 



Eo(x t ) + E t 



NP 



E (x t )+E^ 



NP 



57T 



Eo(x t ) + E t 



NP 



6ti 
C W (M W 



4C (x t )+4Cr + D (x t ) + D* 



NP 







(4) 
(5) 

(6) 

(7) 

(8) 

(9) 
(10) 
(11) 



1 For more details about the structure of model III and the available experimental constraints, one 



can see Refs. [14|,H,17|. 
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C 9 (M W ) = 2p Uc (x t ) + AC» P + D (x t ) + D 



NP 




4 
9 



sin 2 9 



w 



10B (x t )-4C (x t )+4C£ 



NP 



C g {M, 



w 



(12) 
(13) 



where x t = mf/M^r, the functions B (x), Cq(x), D (x), Eo(x) and E' Q are the familiar Inami- 
Lim functions which describe the contributions from the W-penguin and Box diagrams in 
the SM, and can be found easily, for example, in Ref. [I3|. The functions Cq P , Dq P , Eq P 
and E'q P in Eq.(|l3|) describe the new physics contributions to Wilson coefficients in model 



III [0 



n NP 
°0 



A' f 



16 



lit 



l-y t (1 - y t ) 2 



• |A tt | 



NP 



D 



pNP 
E 



-~H(y t )\\ tt \ 2 



with 



H{y) 
i{y) 
J(y) 

K(y) 



NP 
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J(y t )\\tt\ 2 - K(y t )\\ tt \ bb \e 



w 



38y - 79y 2 + 47y 3 Ay - 6y 2 + 3y 4 

+ — rr— h%] 



72(1 -y) 3 
16y - 29y 2 + 7y 3 

36(1 -y) 3 
2y + 5y 2 - y 3 

4(1 — y) 3 

-3y + y 2 



+ 



12(1 -2/) 4 
2y - 3y 2 



+ 



6(1 -y) 4 
3y 2 



log[y] 



2(1 -y)* 



log[y] 



4(1 2(1 -y) 3 



log[y] 



(14) 

(15) 
(16) 
(17) 

(18) 
(19) 
(20) 
(21) 



where x t = m 2 /M^, y t = m 2 /Mjj + , and the small terms proportional to m 2 /m 2 have been 
neglected. 

Since the heavy charged Higgs bosons appeared in model III have been integrated out 
at the scale M w , the QCD running of the Wilson coefficients Ci(M w ) down to the scale 
fi = 0(m b ) after including the NP contributions will be the same as in the SM. In the NDR 
scheme, by using the input parameters as given in Eqs.@ and ( ^7)) and setting \i = 2.5 
GeV, we find that: 



G x = 1.1245, C 2 = -0.2662, C 3 = 0.0186, C A = -0.0458, 
C 5 = 0.0113, C 6 = -0.0587, C 7 = 0.0006, C 8 = 0.0007, 
C 9 = -0.0096, C10 = 0.0026, Cf f = 0.3364 



(22) 



where C e J f = C 8G + C< 
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In this letter, the generalized factorization ansatz as being used in Ref. || will be em- 
ployed. For the studied seven B meson decay modes, we use the decay amplitudes as given 
in Ref. [0] without further discussion about details. We focus on estimating the new physics 
effects on these seven measured decay modes. In the NDR scheme and for SU(3)c, the 
effective Wilson coefficients can be written as [B| 



a 



eff 




T i m b 

+ lv lo S — 



Cj + (C t + C p + C g ) + ^ B>C e , (23) 

ij 



where A[ = (0, 0, -1, 3, -1, 3, 0, 0, 0, 0) T , B[ = (0, 0, 0, 0, 0, 0, 1, 0, 1, 0) T , the matrices ry and 
7v contain the process-independent contributions from the vertex diagrams. The matrix 7y 
and ry have been given explicitly, for example, in Eq.(2.17) and (2.18) of Ref. || 0. The 
function C t , C p , and C g describe the contributions arising from the penguin diagrams of the 
current- current Qi^, the QCD operators Q3-Q6, and the tree-level diagram of the magnetic 
dipole operator Qsg, respectively. The explicit expressions of the functions C t , C p , and C g 



can be found for example in Refs. |8|JlI]| . We here also follow the procedure of Ref. [ jl| to 
include the contribution of magnetic gluon penguin. 

In the generalized factorization ansatz, the effective Wilson coefficients Cf will appear 
in the decay amplitudes in the combinations, 

a^^C^ + ^j, a 2l = C e J f + ^}, (< = 1,...,5) (24) 

iVc iVc 

where the effective number of colors N%" is treated as a free parameter varying in the range 
of 2 < N^* < 00, in order to model the non-factorizable contribution to the hadronic 
matrix elements. It is evident that the reliability of generalized factorization approach has 
been improved since the effective Wilson coefficients C\ appeared in Eq . (|24]) are now 
gauge invariant and infrared safe [I9|. Although iVj" can in principle vary from channel 



to channel, but in the energetic two-body hadronic B meson decays, it is expected to be 
process insensitive as supported by the data f|. 

In the B rest frame, the branching ratios BiB — > PP) can be written as 

B(B -> XY) = ^~^\M(B - XF)| 2 , (25) 

1 tot oTTlVl B 

where T tot (B-) = 3.982 x 1(T 13 GeV and T tot (B° d ) = 4.252 x KT 13 GeV obtained by using 
t(B~) = 1.653ps and t(B%) = 1.548ps PP| , ps is the four-momentum of the B meson, 
M B = 5.279 GeV is the mass of B u or B d meson, and 



is the magnitude of momentum of particle X and Y in the B rest frame. 
In the numerical calculations the following input parameters will be used: 



2 The correct value of the element {rNDR)m and {tndr)%% should be 17 instead of 1 as pointed in 
Ref. |0|. 



5 



The coupling constants, gauge boson masses, light meson masses, 
unit of GeV ) 0,|O| 



(all masses in 



a. 



0.118, sm 2 6 w 



= 1/128, a s (M z 
M z = 91.188, M w = 80. 12. ni Jfl 
m w o = 0.135, m, 



"'/ 



0.547, m v > 



0.23, Gp = 1.16639 x 10" 
= 5.279, m n ± = 0.140, 
0.958, m K ± = 0.494, m K o = 0.498. 



\GeVy 



(27) 



The elements of CKM matrix in the Wolfenstein parametrization: A = 0.81, A = 
0.2205, p = 0.12 and w = 0.34 (which corresponds to 7 = 71° and (3 = 26°), and the 
uncertainty of 5w = ±0.08 will be considered. 

We firstly treat the internal quark masses in the loops as constituent masses, 

m b = 4.88GW, m c = l.hGeV, m s = 0.5GeV, m u = m d = 0.2GeV. (28) 

Secondly, we use the current quark masses for m 8 (i = u,d, s,c,b) which appear through 
the equation of motion when working out the hadronic matrix elements. For p = 
2.5GeV, one finds 

m h = 4.88GeV, m c = 1.5GeV,m s = 0.122GeV, m d = 7.6MeV, m u = 4.2MeV. (29) 

For the mass of heavy top quark we also use mt = wT[(mt) = 168Ge^. 
The decay constants of light mesons (in the units of MeV) are 

fn = 133, f K = 158, /« = f d q = 78, ft = f} = 68, 

ft = -0.9, ft = -0.23, ft = -113, ft = 141. (30) 

where f u ,,-. and f s (/) have been defined in the two- angle-mixing formalism with 6$ = 
-9.1° and 9 8 = -22.2° |ET|. 



The form factors at the zero momentum transfer are 

F B7r (0) = 0.33, F^ K (0) = 0.38, F B,? (0) = 0.145, F W (0) = 0.135 (31) 



in the BSW model ITBl, and 



F B "(0) = 0.36, F£ K (0) = 0.41, Fo V {0) = 0.16, F W (0) = 0.145, (32) 

in the LQQSR approach [[?[]. Here the relation between Fq V (0) and ^^(0) as defined 
in Eq.(Al2) in Ref. has been used. The momentum dependence of F (k 2 ) as defined 
in Ref. |16| is F (k 2 ) = F (0)/(1 — k 2 /m 2 (0 + )). The pole masses being used to evaluate 
the /c 2 -dependence of form factors are m(0 + ) = 5.73 GeV for ub and db currents, and 
m(0 + ) = 5.89 GeV for sb currents. 
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For the seven studied B meson decay modes, currently available measurements from 
CLEO, BaBar and Belle Collaboration Pfl follows: 

BiB -> + - [ (43±}: ' ± °' 5) X 10_6 [CLE0] ' (33) 

7T ) - j (9 3+2.8+1.2) x 10 _ 6 [BaBar]; (33) 

k(b - ir^ - J ( 1L6 -^-i:3) x 10-6 [ CLE °]' ^ 

7T ) - j (lg g± 5.5 ± 2 3) x 1Q _ fl [BeUe]j (34) 

f (17.2l|l ±1.2) x 1(T 6 [CLEO], 
B(B -> K+tT) = i (12.5j£g ij"? ± 2.3) x 1(T 6 [BaBar], (35) 
( (17.4j£J ± 3.4) x 1(T 6 [BELLE], 

B{B -> K°7T + ) = (18.21^0 ± 1-6) x 1(T 6 [CLEO], (36) 

' — | (21±f;l±|;|) x 10- 6 [BELLE], {S7) 

BiB - JTV1 - \ {8 ° + - l ° ± ?) X 10 ^ [CLE ° ] ' (38) 
77 J - \ (62 ± 18 ± 8) x 10- 6 [BaBar], W 

B(B -> J^V) = (89±i| ± 9) x 1(T 6 [CLEO]. (39) 

The measurements of CLEO, BaBar and BELLE Collaboration are consistent with each 
other within errors. 

In Table Q, we present the theoretical predictions of the branching ratios for the seven 
B decay modes in the framework of the SM and model III by using the form factors from 
Bauer, Stech and Wirbel (BSW) model [TO and Lattice QCD/QCD sum rule (LQQSR) 



model p2| , as listed in the first and second entries respectively. The last column shows the 
CLEO data of B — > K°tt + , K°rj' decays, and the average of CLEO, BaBar and/or BELLE 
measurements for other five decay modes. Theoretical predictions are made by using the 
central values of input parameters as given in Eqs.(j^[2~7|-j32|), and assuming M H + = 200GeV 
and N%f* = 2, 3, oo in the generalized factorization approach. The branching ratios collected 
in Table | are the averages of the branching ratios of B and anti-5 decays. The ratio 5B 
describes the new physics correction on the decay ratio and is defined as 



B{B -> XY) IH - B{B -> XY) SM 
BiB -> XYY 



SB(B - XY) = „/ D ( 4 °> 



The last column in shows the CLEO data of B — > K°tt + , K°rj' decays, and the average of 
CLEO, BaBar and/or BELLE measurements for remaining four decay modes. 
From Table |, we find that 

• The SM prediction of B° — ► tt + tt~ decay is clearly much larger than the CLEO mea- 
surement, but agree with BaBar measurement. The new physics contribution to this 
tree-dominated decay mode, however, is negligibly small. 

• For four B — > Kir decays, the SM predictions agree with experimental measurements 
within errors. In model III, the new physics enhancements are large in magnitude: 
~ 50% w.r.t the SM predictions. The model III predictions are generally larger than 
the data for B — > K + tt, K°tc + decays but still agree with the data within 2a errors 
since both the theoretical and experimental errors are still large now. 
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• For B — > Krf decays, the new physics enhancements in model III are large in magni- 
tude: ~ 60% w.r.t the SM predictions. Such enhancement can make the theoretical 
predictions become consistent with the CLEO/BaBar data within one standard devia- 
tion, as illustrated in Fig.p] where the dot-dashed and solid curve shows the theoretical 
prediction in the model III for N^** = 3, oo, respectively. 

• Since the form factors in LQQSR approach are larger than those in the BSW model, 
the theoretical predictions in the LQQSR approach are generally larger than those in 
the BSW model by ~ 15%. 

Because the branching ratios of the studied B decay modes strongly depend on the values 
of involved form factors i^ 7r (0), F^ K (0) and F Bv (0), any information about these form 
factors from data will help us to refine the theoretical predictions. Since the B° — > 7r + 7r~ 
decay is a tree-dominated decay mode and the possible new physics effect is also negligibly 
small, the experimental measures of this decay lead to a stringent constraint on the form 
factor F B7r (0). If we take the average of CLEO and BaBar measurements, 

B(B° d -> 7T+7T-) = (5.5 ± 1.5) x 1(T 6 , (41) 

as the experimental result, then the constraint on F^ n (0) will be 

F B7r (0) = 0.25 ±0.03 (42) 

by setting A = 0.2205, A = 0.81, p = 0.12, 77 = 0.34, N e Jf = 3, and by neglecting FSI also. 

Theoretically, small form factor F^ n (0) will lead to small predicted branching ratios of 
B — > Kir and B — > Krj' decays. First, the relation between Fq V (0) and _F B7r (0) as given in 
Ref. is 

F Bri '(0) = F B ^(0)(sm6 8 /VQ + cosOq/Vs) (43) 

with O = -9, 1° and 9 8 = -22.2°. A small F B7r (0) leads to a small Fq V ' and in turn small 
branching ratios of B — > Krj' decays. Second, F^ K (0) cannot deviate too much from F^ n (0), 
otherwise the SU (3)-symmetry relation Fq 71 »s Fq K will be badly broken. In Table H, we 
show the branching ratios of seven studied decay modes obtained by using F Bjr (0) = 0.25 
instead of 0.33 while keep all other input parameters remain the same as being used in Table 

ffl- 

Contrary to the case of using F^ n (0) = 0.33 in the BSW model, where the inclusion 
of new physics contributions in the model III will degenerate the agreement between the 
theoretical predictions and the data for first three B — ► decay modes, the inclusion of 
~ 50% new physics enhancements to B — > Kn decays for the case of using F B7r (0) = 0.25 
does improve the agreement between the theory and the data, as illustrated in Figs.(|]|l) 
where the short-dashed and solid curve shows the predictions in the SM and model III for 
the case of using F B7r (0) = 0.25. The horizontal band between two dots lines corresponds 
to the (averaged ) data with 2a errors. 

For the decay B -> K n° we find B{B -> K°n°) = (4.3 ± 2.1) x 10" 6 x (F Bw (0)/0.25) 2 
in the SM, which is almost four times smaller than the central value of the averaged data: 
BiB — > K°7i°) = (16.6 ± 5.3) x 10~ 6 . The sixty percent new physics enhancement will be 
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helpful to increase the theoretical prediction, but is still not large enough to cover the gap, 
as illustrated in Fig. (3b) where the lower short-dashed and solid curves show the theoretical 
prediction in the SM and model III with .£^(0) = 0.25. This problem will become clear 
when more precise data from B factories are available. 

For F B7r (0) = 0.25, the SM predictions for branching ratios of B — * Krf decays are in 
the range of (18 — 40) x 10~ 6 as shown in Table |I| and clearly much smaller than the data. 
We know that the Kr( decay rates can be enhanced, for example, through (i) constructive 
interference in gluonic penguin diagrams, which is qualitatively OK but numerical problems 
remain; (ii) the small running mass m s at the scale m& (iii) larger form factor Fq 11 (0) due 
to the smaller rj — rj' mixing angle —15.4° rather than « —20°; (iv) contribution from the 
intrinsic charm content of rj' | 23fl . However, as pointed out in Ref. ||18|| , the above mentioned 



enhancement is partially washed out by the anomaly effects in the matrix element of pseudo- 
scalar densities, an effect overlooked before. As a consequence, the net enhancement may be 
not large enough. For a smaller F^ n (0) = 0.25 ±0.03 preferred by the data, the discrepancy 
between the data and the SM predictions for B — > Krj' decays becomes larger. 

In the model III, however, the new gluonic and electroweak penguin diagrams contribute 
to the B —>■ K + rj' and K°r]' decays through constructive interference with their SM counter- 
parts and consequently provide the large enhancements, ~ 60% w.r.t. the SM predictions, 
to make the theoretical predictions become consistent with the data even for i^^O) = 0.25 
instead of 0.33 as shown in Table || and Fig.f|. 

In Fig.|], we plot the mass-dependence of B(B + — > K + i]') and B(B° — > K°r]') in the 
SM and model III by using F S7r (0) = 0.25 instead of 0.33 (while all other input parameters 
are the same as in Fig.|l|). The short-dashed line in Fig.[| shows the SM predictions with 
N%" = 3. The dot-dashed and solid curve refer to the branching ratios in the model III for 
N^ff = 3 and oo, respectively. The upper dots band corresponds to the data with 2a errors: 
B(B -> K+rf) = (75 ± 20) x 10~ 6 and B(B -»■ K°t]') = (89l$) x 10~ 6 . 

By comparing the curves in Fig.|I] and Fig.[|, it is easy to see that (a) the gap between 
the SM predictions of B — > Krj' decay rates and the data is enlarged by using F c P 7r (0) = 0.25 
instead of 0.33; (b) the new physics enhancement therefore becomes essential for the theo- 
retical predictions to be consistent with CLEO/BaBar result within one standard deviation. 

We know that the calculation of charmless hadronic B meson decay rates suffers from 
many theoretical uncertainties ]7|,|8| . Most of them have been considered in our calculation. If 
we consider effects induced by the uncertainties of major input parameters 77 = 0.34 ± 0.08, 
k 2 = m 2 b /2 ± 2 GeV 2 , F B7r (0) = 0.25 ± 0.03, F<f*(0) = 0.30 ± 0.05, < l/N e Jf < 0.5, 
M H + = 200 ± 100 GeV and m s = 0.1 - 0.122 GeV, we find numerically that 

B (^^)^(^w)={g:?^;°::;:™ delm . (44, 

The SM prediction of B — > Krj' is much smaller than the data, while the model III prediction 
can be consistent with the data within one standard deviation. This is a simple and plausible 
new physics interpretation for the observed rj'K puzzle. 



'However, a rather small m s is not consistent with recent lattice calculations. 
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As simple illustrations we show explicitly the dependence of the branching ratios BiB — > 
Ki]') on the form factor F^ K (0) and the running quark mass m s in Figs. (5, 6). In Fig|| 
we plot the Ff K (0) dependence of the ratios B{B -> K + i) and B(B -> K°r}') in the SM 
and model III, assuming M H+ = 200 GeV, F B7r (0) = 0.25 and Ff K (0) = 0.25 - 0.40. The 
short-dashed line shows the SM predictions with iVjf * = 3. The dot-dashed and solid curve 
refer to the branching ratios in model III for N^** = 3 and oo, respectively. The dots band 
corresponds to the (averaged) data with 2a errors. The theoretical predictions also show a 
strong dependence upon F<f *(0): B(B -> K+rf) = 13.8 x 10" 6 and 26.3 x 10~ 6 in the SM 
for N e Jf = 3 and F<f *(0) = 0.25 and 0.40, respectively. 

In Fig.||, we plot the m s -dependence of the ratios BiB — > K + i]') and BiB — > K°rf) in 
the SM and model III with F B7r (0) = 0.25, F SK (0) = 0.33 and M H+ = 200 GeV. The 
short-dashed line shows the SM predictions with N^f = 3. The dot-dashed and solid curve 
refer to the branching ratios in model III for N%** = 3 and oo, respectively. The dots band 
corresponds to the (averaged) data with 2a errors. The theoretical predictions show a very 
strong dependence upon the mass m s : BiB — > K + i]') = 44.8 x 10~ 6 and 18.6 x 10~ 6 in the 
SM for N^ ff = 3 and m s = 0.08 GeV and 0.15 GeV, respectively. 

In short, we here studied the new physics contributions to the seven observed B — ► PP 
decay modes, and made an effort to find a new physics interpretation for the so-called 
rj'K puzzle of B meson decays by employing the effective Hamiltonian with generalized 
factorization. Within the considered parameter space we found that: 

• The new physics enhancement is negligibly small to tree-dominated B — ► 7r + 7r~ decay, 
but can be significant to the penguin-dominated B — > Kir and B — ► Krj' decay modes, 
~ (40 - 70)% w.r.t the SM predictions. 

• The CLEO/BaBar measurement of B — > 7r + 7r~ decay prefers a small F^ n (0): 
Ff*(0) = 0.25 ± 0.03 instead of 0.33 or 0.36 in the BSW and LQQSR form fac- 
tors. A smaller i^ 7r (0) will leads to smaller predictions for other six B — > Kit and 
B — > Kr]' decay modes studied here. The new physics enhancements to B — > Kit de- 
cays are helpful to improve the agreement between the data and theoretical predictions 
for these decays. 

• The new physics enhancements can boost the theoretical predictions of the branching 
ratios BiB — > K + i]') and BiB — > K°r)') to be consistent with the data within one 
standard deviation. This is a simple and plausible new physics interpretation for the 
observed rj'K puzzle. 
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TABLES 



TABLE I. Branching ratios (in units of 10 6 ) of seven studied B decay modes in the SM and 
model III by using the BSW ( the first entries) and LQQSR (the second entries) form factors, with 
k 2 = ml/2, A = 0.81, A = 0.2205, p = 0.12, 77 = 0.34, 6 = 0°, N^f = 2, 3, 00 and M H+ = 200 
GeV. The last column shows the (averaged) data. 
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TABLE II. Branching ratios (in units of 10 6 ) of seven studied B decay modes in the SM and 
Model III by using the BSW form factors with F B7r (0) = 0.25 instead of F^ n (fi) = 0.33, and 
assuming k 2 = m 2 b /2, A = 0.81, A = 0.2205, p = 0.12, 77 = 0.34, iV c e// = 2, 3, 00, 9 = 0° and 
M H + = 200 GeV. The last column shows the (averaged) data. 
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B°^tt + tt- 5.18 5.89 7.42 5.31 6.03 7.60 2.5 2.5 2.4 5.5 ± 1.5 

B+ K+ir° 7.44 8.36 10.4 10.6 12.0 15.0 43 43 44 13.3 ± 2.9 

B° -► K + 7r~ 10.1 11.2 13.6 15.3 17.0 20.8 51 51 53 15.9 ± 2.2 

B + -► K°tt + 11.5 13.4 17.6 17.2 19.9 26.0 50 49 48 18.2±|;§ 

B° -► K°tt° 3.79 4.29 5.46 6.04 6.85 8.70 60 60 60 16.6 ± 5.3 

B+ K+rf 19.1 24.4 36.9 32.5 40.4 58.9 70 66 60 75 ± 10 

B° -► K°7]' 18.3 23.7 36.6 30.9 38.9 57.6 69 64 58 89l^ 
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FIGURES 
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FIG. 1. Plots of branching ratios of decays B + — ► K + rj' (la) and B° — ► i^T ?/ (lb) versus 
mass M H + in the SM and model III with F B7r (0) = 0.33. The short-dashed line shows the 
SM predictions with N%ff = 3. The dot-dashed and solid curve refers to the branching 
ratios in the model III for N%ff = 3 and oo, respectively. The dots band corresponds to the 
CLEO/BaBar data with 2a errors. 
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FIG. 2. Plots of branching ratios of decays B — > K + 7t° (2a) and K + ir (2b) versus 
l/Np* in the SM and model III assuming Mjj+ = 200 GeV. The short-dashed and solid 
curve show the predictions in the SM and model III using F ( f 7r (0) = 0.25. The dots band 
corresponds to the (averaged) data with 2a errors. 
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FIG. 4. Same as Fig.@ but for Ff*(fl) = 0.25. 
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FIG. 5. Plots of branching ratios of decays B + — ► K + rj' (5a) and B — > .ftT ?/ (5b) versus 
form factor F(f ^ (0) in the SM and model III. The short-dashed line shows the SM predictions 
with N%ff = 3. The dot-dashed and solid curve refers to the branching ratios in the model 
III for N^ff = 3 and oo, respectively. The dots band corresponds to the (averaged) data with 
2cr errors. 
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FIG. 6. Plots of branching ratios of decays B + — ► K + rj' (6a) and B — > K°r]' (6b) versus 
mass m s in the SM and model III with F^ n = 0.25, F^ K = 0.33 and M H + = 200 GeV. The 
short-dashed line shows the SM predictions with N%f f = 3. The dot-dashed and solid curve 
refers to the branching ratios in the model III for N%Jf = 3 and oo, respectively. The dots 
band corresponds to the (averaged) data with 2a errors. 
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